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cosine -squared  spreading  function  is  used  to  model  the  sea  state  environment.  An 
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1.  INTRODUCTION 


The  ocean  waveform  is  an  undulating  process  that  transmits  energy  thousands  of 
miles  both  efficiently  and  effectively.  A  twelve-second  periodic  wave  group  will  take 
forty- eight  hours  to  cross  1,000  miles  of  deep,  open  ocean  with  little  loss  of  energy. 
When  its  power  is  unleashed  on  a  rigid  body,  it  becomes  one  of  the  greatest  fears  of  all 
mariners  and  maritime  insurance  companies  alike.  It  is  the  effort  of  the  naval  architect  to 
design  ships  that  can  surpass  a  broad  spectrum  of  waveforms  and  the  energy  associated 
with  them. 

A  ship  advancing  at  a  relatively  constant  forward  speed  with  an  arbitrary  heading 
in  a  train  of  regular  wa’ves  will  move  in  six  degrees- of- freedom  (three  translational: 
surge,  heave,  sway;  three  rotational:  pitch,  roll,  yaw).  In  order  to  be  able  to  compute  all 
responses  of  the  vessels  to  random  seaways,  it  is  necessary  to  deal  with  the  complete 
motions  of  the  vessel  in  all  six  degrees- of- freedom,  with  special  attention  to  the  important 
couplings  of  heave  and  pitch,  sway  and  yaw.  For  such  a  rigid  body,  six  nonlinear 
equations  of  motion,  with  six  unknowns,  must  be  set  up  and  solved  simultaneously.  In 
the  case  of  a  vessel  of  port/starboard  symmetry  and  restrained  about  this  axis  of 
symmetry,  the  six  non-linear  equations  can  be  reduced  to  two  sets  of  three  linear 
equations. 

Specifically,  it  is  the  effort  of  this  thesis  to  present  a  linear  theory  of  ship  motions. 
While  most  vessel  responses  are  nonlinear  to  some  extent,  however,  where  nonlinearities 
are  small  a  linear  theory  will  yield  reasonable  predictions.  Experimental  and  theoretical 
investigations  have  repeatedly  shown  that  a  linear  theory  analysis  gives  excellent 
predictions  over  a  wide  variety  of  sea  states  and  vessel  platforms  (Beck). 
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A.  MOTIVATION 

Since  the  advent  of  steam,  the  level  of  engineering  necessary  to  keep  ocean-going 
vessels  afloat  and  righting  has  had  to  become  more  and  more  complex.  The  economic 
demand  of  trading  nations,  and  the  warship  designs  of  international  navies  continually 
pushes  the  complexity  level  to  yet  another  extreme.  Since  the  rapid  developments  in 
hydrodynamic  theory  of  the  1950s,  the  advances  in  oceanography  and  computer 
technology  have  enveloped  the  marine  engineers  task  in  pursuing  advanced  designs  in 
high-speed  hull  forms,  and  unusual  platforms  such  as  SWATH  vessels. 

Close  proximity  towing  of  surface  ships  is  a  matter  of  interest  to  the  U.S.  Navy 
and  the  Office  of  Naval  Research.  Several  possible  applications  have  been  suggested 
including  the  SLICE/KAIMALINO  connection  (Nash)  and  the  SEA  EANCE 
configuration  (TSSE).  Of  primary  concern  in  this  dissertation  is  the  evaluation  of  motions 
in  six  degrees-of- freedom  of  the  two  ships  when  in  close  proximity  towing  operations. 
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II.  WAVE  SPECTRA 


A.  EARLY  WAVE  THEORY 

Franz  Gerstner  of  Czechoslovakia  officially  documented  the  first,  rather  primitive 
wave  theory  in  1802.  He  described  how  water  particles  move  in  oscillatory  motions 
within  the  wave.  He  further  described  how  the  water  particles  in  the  crest  move  in  the 
direction  of  wave  advance  and  those  in  the  trough  move  in  the  opposite  direction.  Such 
observations  had  been  made  centuries  earlier;  it  was  his  discovery,  however,  that 
described  the  circular  path  of  the  water  particles  to  have  a  diameter  equal  to  the  height  of 
the  passing  wave  (path  diameter  proportional  to  water  depth).  He  essentially  described  a 
trochoidal  wave  (Figure  1).  This  discovery  bridged  the  gap  between  the  casual  seafarer’s 
observation  to  that  of  mathematical  derivation  and  hence  the  theoretical  beginning  of 
hydrodynamics. 

Following  the  theoretical 
came  the  experimental.  The  first 
experimentalists  developed 

controllable  wave  tanks  with  the 
capability  to  reproduce  waveforms 
over  and  over,  and  whereby  careful 
measurements  could  be  taken. 

Then  in  the  1950’s  a  new  breed  of 
experiment  was  developed 
harnessing  modern  computations 
and  statistical  applications;  and, 
hence,  the  numerical  wave  tank  was 
bom. 
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Figure  1.  Trochoidal  Wave  (After  Beck) 
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B.  WAVE  SPECTRA  ANALYSIS 


The  understanding  of  vessel  response  at  sea,  and  the  ability  to  predict  their 
behavior  begins  with  the  study  of  the  nature  of  ocean  waves  in  which  the  vessels  will 
operate.  Consequently,  in  order  to  understand  the  behavior  of  ocean  waves,  one  must 
first  understand  their  origin.  Waves  are  a  natural  phenomenon  that  occurs  between  the 
interface  of  any  two  fluids  of  different  density.  While  any  kind  of  disturbance  in  the 
ocean  will  generate  waves  of  sparing  severity,  there  exist  four  primary  sources:  wind, 
earthquakes,  landslides,  and  the  gravitational  attraction  of  the  moon  and  sun.  Of  the  four, 
the  greatest  of  importance  is  the  wind- generated  waves.  The  size  and  variety  of  these 
wind  waves  depends  on  many  factors,  namely,  the  velocity  of  the  wind,  the  distance  it 
blows  across  the  water,  the  time  duration  that  it  blows,  and  the  depth  of  water  in  which 
the  wind  trans verses  (Figure  2).  The  character  and  signature  of  the  wave  is  quickly  lost 
as  the  wind  energy  is  transferred  to  tidal  energy  in  this  highly  stochastic  ocean  media. 
Once  the  winds  have  transferred  their  energy  to  the  ocean  surface,  the  waves  begin  their 
natural  procession  from  the  generation  area.  The  original  wind  waves  eventually  decay, 
evidenced  by  the  crests  as  they  become  lower,  and  more  rounded  and  symmetrical. 
These  crests  then  begin  to  move  in  groups  of  similar  period  and  wave  height,  and  their 
form  becomes  more  and  more  sinusoidal.  This  origination  theory  becomes  now  the  basis 
for  our  analytic  discussion:  sinusoidal  wave  forms  moving  as  steady- state,  constant 
parameter  wave  “trains”  under  the  control  of  gravity  and  inertia;  and  their  subsequent 
influence  on  sea- going  vessels  in- tow. 


4 


0  0.5  1  1.5  2  2  5 

Wave  Frequency,  rad/sec 

Figure  3.  Sea  State  Four  (Varying  Period)  (From:  McCreight) 


Waves  are  classified  according  to  their  period,  which  can  range  from  a  fraction  of 
a  second  to  more  than  ten  thousand  seconds.  They  can  range  in  size  from  a  ripple  to 
great  storm  waves  and  tides  whose  wavelength  can  reach  half  the  circumference  of  the 
earth  (Bascom).  Essentially  we  are  studying  a  medium  whereby  the  wave  spectrum 
ranges  from  waves  so  small  they  can  barely  be  seen  to  waves  so  long  that  they  are  not 
even  noticed.  For  these  extreme  conditions  our  vessel  design  is  safe.  It  is  a  narrow  band 
of  wave  spectrum  that  can  cause  greatest  damage  to  seagoing  vessels  and  is  of  utmost 
concern  to  the  naval  architect  and  our  research. 

The  energy  in  the  ocean  is  distributed  among  several  distinct  groups  of  waves 
each  with  a  defined  range  of  periods  (Figure  3).  All  waves  can  be  classified  as  gravity 
waves  since  once  they  are  started;  gravity  is  the  driving  force  that  keeps  these  sea  crests 
moving  in  a  never-ending  attempt  to  restore  the  original  flat  sea  surface.  The  potential 
energy  of  each  wave  is  due  in  part  to  the  waves  profile  which  consists  of  two  sides  of  the 
equilibrium  position:  the  crest  that  rises  above  datum  sea  line,  and  the  trough  that 
extends  below.  A  wave  group  moves  perpetually  attempting  to  overtake  the  trough  ahead 
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and  restore  equilibrium,  suffering  of  course  from  the  wave  ahead  with  the  same  intent. 
For  swell,  the  wave’s  potential  and  kinetic  energies  are  equal.  The  kinetic  energy 
possessed  by  a  wave  is  due  to  the  imparted  motion  of  the  water  particles,  while  the 
remainder  of  the  waves  energy  is  the  potential  energy  due  to  the  elevatbn  of  the  center  of 
gravity  of  the  mass  of  water  in  the  crest  above  mean  sea  level  (Figure  4). 


Another  classification  of  waves,  although  when  fully  developed  act  like  gravity 
waves,  are  called  capillary  waves  (short  crested  seas).  They  are  a  special  waveform 
whose  driving  force  is  not  gravity,  but  rather  the  surface  tension  at  the  water/wind 
interface.  The  capillary  force  is  stronger  than  gravity,  and  given  the  open  ocean 
condition,  thrives  in  producing  the  most  abundant  kind  of  wind- generated  wave. 
Capillary  waves  give  rise  to  the  development  of  ripples  (incident  wind  drag),  which  are  at 
the  lower  end  of  the  wave  spectrum,  but  none-the-less  lead  to  the  growth  and  overall 
combined  effect  of  fully  developed  seas  (wave  energy  absorption  proportional  to  wind 
velocity). 

Open  ocean  waves  do  not  have  the  normal  and  precise  properties  as  those 
generated  in  a  wave  generator.  The  height  of  the  crests  and  the  depths  of  the  troughs  are 
highly  irregular.  The  theoretical  terminology  of  a  discernable  wave  period  and  wave 
velocity  (celerity)  is  lost  when  taken  to  the  real  and  observable,  highly  stochastic  and 
irregular  environmental  oceanic  swirl.  For  these  very  reasons  is  why  the  physics  of 
ocean  waves  must  rely  on  statistical  methods  in  order  to  describe  the  properties  of  waves 
(Appendix  B). 

Wave  heights,  while  apparently  just  as  irregular,  can  be  referenced  to  a  datum  line 
also  referred  to  as  mean  sea  level.  Using  mean  sea  level  as  a  baseline,  a  sea  state  is 
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simply  the  result  of  superimposing  a  number  of  sinusoidal  wave  trains  one  on  top  of  the 
other  (Figure  5). 

Each  layer  represents  a  series  of  regular  sine  waves,  each  having  their  own 
distinct  wave  height,  wavelength,  and  direction;  and  coincides  with  the  theoretical 
notions  exactly.  The  summation  of  the  wave  crests  and  troughs  of  varying  amplitude 
tend  to  cancel  each  other  out,  and  conversely,  amplify  the  local  maximum  or  minimum 
depending  on  coincidence  of  discrete  wave  points  throughout  the  ocean  surface.  Again, 
as  these  waves  transverse  each  other,  the  cancellation  and  amplification  will  subside 
quickly  as  these  respective  wave  groups  transverse  their  own  direction.  Accordingly,  the 
greater  the  number  of  wave  layers  the  more  random  the  sea  surface. 

The  seemingly  random  nature  of  ocean  waves  can  be  better  characterized  by  its 
energy  spectra.  This  methodology  assigns  a  value  to  the  square  of  the  wave  height  for 
each  frequency  and  direction.  Then,  using  statistical  methods,  the  wave  spectrum  can  be 
analyzed  to  reveal  how  the  wave  energy  is  distributed  among  various  wave  periods. 


Figure  5.  Wave  Fourier  Composition  (From:  Ochi) 
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The  essential  feature  of  these  mathematical  techniques/models  is  the  concept  of  a 
spectrum  defining  the  distribution  of  energy  among  the  different  hypothetical  regular 
wave  components  having  various  frequencies  (wave  lengths)  and  directions.  It  has  been 
formulated  that  the  irregular  motions  of  a  ship  in  a  seaway  can  be  better  described  as  the 
linear  superposition  of  the  responses  of  the  same  ship  to  all  the  wave  components  of  the 
seaway.  For  our  purposes,  the  linear  theory  assumption,  as  applied  to  the  SWATH 
vessels,  implies  the  vessels’  response  is  described  by  the  same  statistical  properties  of  a 
random  (stochastic)  surface  wave  environment.  Given  these  assumptions,  these  surface 
waves  can  be  characterized  as  having  a  Gaussian,  or  normal,  probability  distribution 
about  some  mean  (zero  amplitude,  calm  seas),  under  short-term  statistically  stationary 
conditions  (Figure  4).  This  statistical  technique  will  greatly  simplify  the  numerical 
involvement  of  our  study  with  regards  to  the  application  of  statistics,  probability  theory 
and  Fourier  analysis  techniques. 

During  the  energy  transfer  between  wind  and  water,  there  are  several  important 
wave  interactions  and  wave -breaking  processes  that  affect  the  dispersion  and  propagation 
of  waves  from  the  storm  area.  Particularly  for  the  case  of  short  crested,  small  amplitude 
waves,  the  principle  of  linear  superposition  applies.  For  example,  if  Xi(xi,yi,ti)  and 
^2(x2,y242)  are  two  wave  systems,  then  Xi(xi,yi,ti)  -i-  )i2(x2,y2,t2)  is  also  a  wave  system. 
Inherent  to  this  assumption,  the  two  systems  can  move  through  each  other  without 
damping  the  other  whereby  all  wave  energy  is  conserved  linearly.  Of  equal  conceptual 
concern  is  wave  celerity,  or  velocity,  which  is  wholly  a  function  of  wavelength  (shorter 
waves  travel  slower  than  longer  waves).  Referring  to  Figure  5,  it  can  be  shown  that  any 
wave  system  can  be  divided  into  a  sum  of  component  regular  waves  of  various 
wavelength,  amplitude  and  direction  using  Fourier  Integral  techniques.  Since  water  is 
considered  to  be  incompressible,  the  average  value  of  vertical  displacement  at  any  instant 
for  a  sinusoidal,  or  regular,  wave  period  is  zero.  While  the  mean  will  give  us  little 
statistical  information,  the  variance,  however,  is  a  positive  value  that  will  directly  relay 
the  severity  of  the  sea.  A  fundamental  theory  of  statistics  states  that  the  variance  of  the 
entire  system  is  determined  by  measuring  the  amount  of  variability  in  the  distribution  of 
the  independent  wave  groups  themselves.  In  general,  the  variance  of  a  continuous 
function  can  be  expressed  by  the  following: 
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And  in  the  case  of  simple  harmonic  sine  waves,  the  variance  of  wave  elevation  of 
a  single  cycle  is  equal  to  one-half  the  square  of  the  wave  amplitude. 

The  slight  errors  involved  in  the  assumed  sinusoidal  wave  shape  of  linear 
harmonic  wave  theory  are  standard  practice,  and  can  be  considered  reasonable 
assumptions  with  little  consequence.  The  errors  only  become  significant  when  dealing 
particularly  in  the  wave  breaking  phenomena,  and  when  the  geometry  of  the  wave  crests 
becomes  unsteady.  The  correct  mathematical  analysis  of  linear  short- crested  irregular 
waves  is  formable  with  little  error  using  short-term  stochastic  models  of  ocean  waves. 
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^  (x,y,  t)-^  cos[k(vcos  p  -h  ysin  P  )-(0t-l-8] 


8=phase  angle 

P=wave  direction  (global  coordinates) 
(D=wave  frequency 
Wave  Number:  k  =  2k  I L^,  -  I g 

y2 

Wave  Length:  =  271  — ^ 

s 


Wave  Period:  T  = 


^  2kL  ^ 


^  J 


Surface  Profile:  Co  =  C  cosk{x-Vj) 


Velocity  Potential:  ^  =-^V^e^  smk(x-Vd) 


Wave  Variance: 


Wave  Celerity: 


y 


T 


'-f 

u 


(Appendix  B) 
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Continuing  with  the  short-term  model,  the  total  wave  system  can  be  shown  as  a 
summation  of  many,  if  not  infinite,  independent  wave  components: 

C(0=LC;cos(-(D/  +  8,)  (2) 

i 


This  system  can  also  be  represented  by  a  point  spectrum,  also  known  as  a 
variance  spectrum,  S((0)  (the  wave  is  observed  only  at  a  point,  and  none  with  respect  to 
wave  direction).  At  any  particular  wave  frequency,  (Di,  the  variance  of  all  the  wave 
components  within  a  small  finite  frequency  band,  5(0,  centered  at  (Di  is  denoted  by: 

(C")^5((0,)5(0  (3) 

following: 

(v)=y2v  (4) 

substituting: 

=  S((oy8(0  (5) 

Integrating  the  RHS  of  the  equation  will  resolve  a  good  approximation  of  the 
wave  system’s  total  energy  (Figure  3). 


£’=|s((o)d(0  (6) 

0 


In  the  derivation  of  energy  spectra  in  short  crested  seas.  Figure  6  depicts  only 


fourteen  wave  groups,  a  multi- directional  short  crested  sea  requires  many  more  (73  was 


icai  E  OF  Tiwe  OR  OtS'AMCE 

Figure  6.  Variance  Spectmm  of  Waves  (From:  Beck) 
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Figure  6,  depicts  several  distinct  sine  curves  of  random  phase.  Numerical 
summing  these  will  result  in  the  irregular  wave  pattern  in  which  we  are  interested.  Figure 
5  demonstrates  this  idea  quite  clearly. 

While  this  seems  a  satisfactory  depiction  of  the  development  of  ocean  wave 
spectra,  and  mathematically  the  variance,  E,  obtained  from  a  point  spectrum  is  a  decent 
measure  of  the  sea  severity,  it  can  be  further  shown  that  this  measurement  is  still 
insufficient.  A  more  complete  characterization  of  the  seaway  must  include  wave 
component  direction,  also  known  as  a  directional  spectrum.  From  the  two-dimensional 
case  to  the  three-dimensional  case  our  generalized  wave  equation  becomes: 

C(x,y,t)  =  IE?  ^  cos[k,(xcos  -I- ysin  P^.)-cqt-i-8,y]  (7) 

i  j 

and 

(c  (t)"}  ^  =  F  F  S((o  ,|3  (8) 

1  J 

Wave  energy  for  a  point  has  an  angular  distribution  as  well  as  a  distribution  over 
a  range  of  frequencies.  This  angular  distribution  of  wave  energy  is  termed  directional 
spreading  and  overall  greatly  increases  the  accuracy  of  our  results  when  designing  for 
short  crested  sea  behavior.  The  directional  spectrum  thus  defines  the  distribution  of 
energy  with  even  greater  resolve  defined  by  the  frequency  and  direction  of  the  variances 
of  the  individual  wave  system  components  (each  having  a  unique  combination  of 
frequency  and  direction,  and  random  phase  angle).  Furthermore,  for  our  incorporated 
spreading  function,  the  directional  spectrum  can  be  subdivided  into  the  following: 

5((0,P)  =  S((0)M(P)  (9) 

where: 

M(p)  =  -cos'p  (10) 

7t 

While  the  frequency  and  direction  are  accounted  for,  what  is  critical  to  our  study 
is  defining  an  idealized  random  seaway.  We  have  already  discussed  the  constructs  of 
wave  spectra  from  recorded  data,  but  how  can  we  reconstruct  such  data  given  certain 
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environmental  parameters?  This  technique  is  particular  to  the  design  process,  where  we 
are  only  interested  in  recreating  wave  energy  spectroscopy.  This  approach  allows  the 
designer  to  modify  such  environmental  parameters  incrementally  to  discern  ship  response 
to  specific  wave  parameters. 

The  Bretschneider  Spectrum  was  selected  and  can  be  presented  as  follows: 

=  (11) 

4  (0 


2n 

T 

m 


(12) 


where: 


Tm  =input  modal  wave  period 

Hs=input  significant  wave  height  defined  as  the  average 
one-third  highest  wave  heights  (H1/.3) 

This  is  a  two-parameter  wave  spectrum  and  is  a  widely  accepted  design  parameter 
for  marine  systems.  Using  this  infamous  two-parameter  spectral  formulation  it  is 
possible  to  generate  a  myriad  of  wave  groups  of  varying  severity  and  size  composition. 
This  spectral  density  function  can  quickly  illustrate  the  magnitude  of  wave  energy.  The 
most  commonly  used  definition  of  sea  severity  is  significant  wave  height,  which  can  be 
deduced  from  the  wave  spectrum  as  being  equal  to  four  times  the  square-root  of  the  area 
under  the  curve  (Figure  3). 

Continuing  with  the  linear  theory  assumption,  the  spectral  density  of  any  response 
can  be  found  by  multiplying  the  incident  wave  spectrum  by  the  square  of  the  response 
amplitude  operator  (RAO)  of  the  desired  response.  In  control  theory,  for  example,  the 
RAO  is  often  called  the  transfer  function  of  the  linear  system  function.  For  any  given 
frequency,  the  RAO  is  the  amplitude  and  phase  of  the  desired  response  to  regular  incident 
waves.  Also  known  as  the  motion  transfer  function,  the  response  amplitude  operator 
(RAO)  maps  the  complex  response  of  a  vessel  to  a  seaway  or  input  spectrum  as  a 
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function  of  frequency.  The  spectrum  of  the  particular  response  that  the  RAO  has  been 
computed  for  is  calculated  from: 

S«(a))  =  |RA6>(a))|"s(a),13)  (23) 

In  this  equation,  5^(co)is  the  response  of  the  vessel  to  the  input  sea  spectrum  for 

a  given  frequency.  Ship  response  then  can  be  calculated  for  a  given  wave  frequency  and 
significant  wave  height.  For  example,  the  complex  absolute  motions  predicted  in  regular 
wave  modeling  (?  §,  ?k)  are  converted  into  RAO’s  for  absolution  motion  as 

=  (24) 

The  response  spectrum  for  absolute  motion  then  becomes: 

‘5«-^.,.(03)  =  |«Mi,,)f>S(to,13)  (25) 

Finally,  all  motion  analysis  speed  polar  plots  are  resolved  in  RMS  values: 

(RM5)"  (26) 

P  (0 
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m.  SHIP  RESPONSE 


The  longitudinal  motions  of  a  symmetrical  ship  with  symmetrical  moorings  in 
regular  waves  can  be  considered  distinctly  from  the  horizontal  degrees  of  motion.  The 
principle  modes  are  the  heave  and  pitch  degrees- of- freedom  (modes  3  and  5, 
respectively).  In  particular,  the  cross  coupling  between  these  two  modes  is  critical  to  the 
success  of  this  model.  To  further  this  discussion,  it  has  also  been  shown  that  the  surge 
degree  of  freedom  (mode  1),  given  the  SWATH  vessels’  architecture,  has  only  a  minor 
contribution  and  can  be  neglected  for  the  present  discussion  (refer:  Conclusions  pp.23) 
For  the  following  derivatbn,  the  SLICE/KAIMALINO  coupling  is  rigid,  has  constant 
forward  velocity,  and  are  resolved  into  the  simplified  head  seas  case  (waves  from  the 
bow;  |i=180°).  It  is  also  assumed  that  the  wave  excitation  forces  are  linear  and  harmonic, 
acting  according  to  the  wave  encounter  frequency  Poppler  Shift).  Using  the  linear 
theory  assumption,  the  responses  of  the  vessel  will  be  directly  proportional  to  wave 
amplitude  and  the  frequency  of  occurrence  corresponds  to  the  frequency  of  encounter. 
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From  Newton’s  second  law  of  motion,  a  translational  displacement  is  defined  as 
the  mass  times  the  acceleration,  therefore,  for  heave  we  have  (center  of  gravity  located  at 
the  waterline): 

Ati-=F3  (27) 

and  for  pitch,  the  rotational  modes  equal  the  mass  moment  of  inertia  times  the  angular 
acceleration: 

(28) 

where  A  equals  displacement,  I55  is  the  mass  moment  of  inertia  about  the  y-axis  and  F3 
and  F5  are  the  total  force  and  moment  acting  on  the  body. 

For  our  discussion  (linear  theory),  the  total  force  and  moment  consist  mainly  of 
hydrostatic  and  hydrodynamic  fluid  forces: 

F3(t)=Fex3(t)+FH3(t)  (29) 

and 

F5(t)=Fex5(t)+FH5(t)  (30) 

where  Fex  is  the  exciting  force  due  to  the  waves  acting  on  the  restrained  ship,  and  Fh  is 
the  radiation  force  due  to  ship  motions  in  an  assumed  calm  sea,  ideal  case.  Again  due  to 
a  linear  theory  assumption,  the  fluid  forces  can  be  divided  as  such. 

F£X3(0  ~  |cos(C0gt +£3)  (31) 

F£X5(^)  ~  ^5)  (32) 

where  IFex3I  refers  to  the  amplitude  of  the  heave  force  and  IFexsI  refers  to  the  amplitude 
of  the  pitching  moment.  Also  called  the  Froude- Krylov,  or  Korvin-Kroukovsky, 
excitations  (NASH),  these  excitations  represent  the  integration  of  the  pressure  field,  over 
the  ship  wetted  surface  area,  which  would  have  existed  in  the  incident  wave  system  if  the 
ship  were  not  present.  The  diffraction  excitation  forces  and  moments,  which  are  also  a 
component  of  IFex3I  and  IFexsI,  are  caused  by  the  diffractbn  or  modification  of  the 
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oncoming  waveforms.  The  respective  epsilon  values  are  the  phase  angles  between  the 
excitation  force  and  sea  waves. 

Pm  =  +  ^33(00)^3  +  A35(a))fi5  +  +  C35(a))ri5  (33) 

Pm  =  [43(«>)n'3  +  553(a))p3  +C53((d)P3  +  45(a) )fi'5  +  555(a))p5  +  Cg5(o))ri5  (34) 

where  Ajk-terms  correspond  to  added  mass  in  phase  with  vertical  accelerations,  the  Bjk- 
terms  to  hydrodynamic  damping  in  phase  with  vertical  velocity,  and  the  C]k-terms  to 
restoring  forces  (buoyancy  effects)  in  phase  with  ship  displacements.  These  coefficients 
were  developed  using  a  strip  theory  (a  fast  and  reliable  numerical  method,  capable  of 
accommodating  a  wide  range  of  hull  forms).  ( A3g(a))rj5  represents  the  force  in  the  heave 
mode  due  to  an  acceleration  in  the  pitch  mode.) 

As  may  be  evident  to  the  reader,  the  coupled  equations  are  similar  in  form  to  that 
of  a  two -degree- of- freedom  mass- spring- damper  system  (added  mass-i-mass,  damping, 
and  restoring  terms  on  the  LHS;  excitation  forces  on  the  RHS).  The  dynamics  of  the  two 
systems,  however,  are  far  from  similar.  For  the  ship  system  all  coefficients  and 
excitation  forces  are  all  functions  of  frequency,  whereas  for  the  simple  mass- spring- 
damper  system  these  values  are  time  independent.  While  solutions  to  the  vessel  coupling 
problem  are  possible  in  the  time  domain,  such  methods  require  heavy  mathematics 
involving  complicated  convolution  integrals.  In  order  to  overcome  this  difficulty  the 
problem  need  only  be  solved  in  the  frequency  domain.  The  exciting  force,  therefore,  can 
now  be  defined  in  the  frequency  domain  as: 

Fe.(t)=F,,e”«‘  (35) 

Defining  the  motion  and  its  derivatives: 

ri(0  ='ne'“'' 

tKO  =  (36) 

m  = 

where  is  the  complex  response  amplitude  containing  both  magnitude  and  phase  of  the 
response. 
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Finally,  linearizing  the  equations  of  motions  (Abkowitz,  1969),  (note  Figure  1): 


A(Ti;  +  z,Ti'5)  =  i"i 
A(Ti2  -  Zjf],  +  =  F, 


A(Ti3-ATl5)  =  ^3 

^44^14  ~  ^46^6  ~  Az/i2  “  ^4 
^55^5  aIz/Ii  ~  “  ^5 


(37) 


^66^6  ^6^4  A-r^'n2  “  7^6 


Rewriting  the  equation  of  motion  in  generalized  form: 


£[-0),'(A,.  +  A,,)  +m,Bi,  +  =  F/  +F,'> 


(38) 
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IV.  RESULTS 


The  following  results  (Table  1)  are  presented  in  speed  polar  plots  (Appendix  A). 
As  discussed  earlier,  waves  on  the  bow,  or  head  seas,  correspond  to  180°  on  the  plots. 
The  input  parameters  as  necessitated  by  the  Bretschneider  Equation  are  modal  period  and 
significant  wave  height,  with  the  additional  input  of  connection  ratio  on  behalf  of  vessel 
coupling  t)  determine  connection  force  and,  ultimately,  modal  excitation.  The  radial 
coordinates  range  from  zero  to  twenty  knots,  and  all  results  show  contours  of  constant 
connection  force  (RMS  values). 

Table  1 .  Numerical  Test  Matrix 


g) 

'CD 

I 


0 

> 

CO 


c 

CO 

o 


c 

g) 

lo 


Modal  Period 

Connection  Ratio 

5 

5 

0.1 

0.5 

1.0 

5 

10 

0.1 

0.5 

1.0 

b. 

15 

0.1 

0.5 

1.0 

10 

8 

0.1 

0.5 

1.0 

10 

12 

0.1 

0.5 

1.0 

10 

16 

0.1 

0.5 

1.0 

20 

8 

0.1 

0.5 

1.0 

20 

12 

0.1 

0.5 

1.0 

20 

16 

0.1 

0.5 

1.0 

30 

12 

0.1 

0.5 

1.0 

3Q 

16 

0.1 

0.5 

1.0 

30 

20 

0.1 

0.5 

1.0 

Somewhat  illusive  to  the  discussion  has  been  from  where  are  such  motions 
observed?  Up  to  this  point  it  has  been  assumed  that  the  wave  spectra  has  been  with 
respect  to  some  fixed  referenced  point.  The  absolute  displacements,  velocities,  and 
accelerations  have  all  been  assumed  to  be  with  respect  to  an  absolute  or  non- 
rotating/transversing  observer.  However,  if  the  velocity  of  the  wave  groups  along  with 
the  velocity  of  the  vessel  is  taken  into  account  the  results  will  change  significantly.  The 
time  record,  or  frequency  of  encounter,  of  a  moving  point  relative  to  a  moving  fluid 
media  and  some  global  fixed  inertia  reference  frame  will  undoubtedly  be  significantly 
affected  by  the  Doppler  Shift  in  the  component  frequencies  of  the  wave  system  itself. 
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I 


Figure  8.  Frequency  of  Encounter  (From:  Ochi) 


The  location  of  the  ship  can  be  given  by: 


=x+Uj 
yo  =  y 


(40) 


Where  the  component  Uot  is  the  difference  between  the  body  fixed  frame  and  the 
global  coordinates.  Substituting  into  the  two-dimensional  wave  equation,  the  expression 
for  the  wave  field  as  observed  from  the  moving  ship  becomes: 

^  ^  cos[/rxcos  |i-i-  /tysin  p  - (co-  kU^cos  p)t  -i-e]  (41) 


From  Figure  8  the  frequency  of  encounter  can  then  be  described  as: 


I 

U„-U„  cosp 


(42) 


W.  =  -^(^w-^oCosp)  (43) 

CO^ 

=  co„-— t/„cosp  (44) 

s 

co^  =  co(l-a)  (44) 
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For  the  simplified  solution: 

a<0:  bow  waves,  (De>0) 

0<a0.5:  stem  waves  a)e<0).  (waves  pass  the  ship  quickly) 
maximum  (Oe  when  a=0.5 

0.5<a<l:  stem  waves  a)e<0)  (waves  pass  the  ship  slowly) 
a=l:  wave  velocity  equal  to  that  of  the  ship. 

a>l:  waves  come  from  the  stern,  but  the  ship  overtakes  the  waves  so 
that  the  waves  appear  to  actually  come  from  ahead. 

The  change  in  frequency  from  that  of  oscillation  to  wave  encounter  is  analogous 
to  the  Doppler  Shift  found  in  sound  waves  and  electro- magnetic  theory.  Considering 
ahead  seas,  the  frequency  of  encounter  will  be  greater  than  that  of  the  absolute  frequency 
and  is  therefore  another  critical  ingredient  to  the  physics  of  vessel  motion  and  overall 
understanding  of  our  results.  The  frequency  of  encounter  and  according  frequency  shifts 
attributed  to  vessel  velocity  has  a  profound  influence  on  ship  motions.  In  particular, 
noting  the  circular,  or  concentric,  motion  amplitudes  with  high  and  low,  and  then  high 
again  motion  amplitude  results.  These  results  are  due  in  part  to  the  speed-polar  plot  data 
representation.  With  zero  velocity  at  the  center,  and  increasing  outwardly  in  integer 
intervals  to  a  maximum  velocity  of  twenty  knots,  the  frequency  shifts  are  rather  evident. 
All  statistical  assumptions  remain  the  same,  as  only  the  point  of  reference  has  changed; 
however,  the  directional  spectrum  must  also  be  slightly  modified  to  account  for  such  said 
frequency  shifts. 

= rr  p  («) 

j 

Realizing  the  spectmm  of  encounter  to  be  the  limit  as  d  P  dco  tends  to  zero: 

m)rP.]  =  >2^^Sa),6p  (46) 

such  that: 
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5(0  „ 


1- 


2(0 1/„  cos  (3 


5(0  (47) 


solving: 

5(0),,  p )  =  5(0) , (3 )/[  1  -  (2(0t/,  /  g  )cos  p  ]  (48) 

The  null  space  on  the  plots  is  a  combination  of  when  a>l,  and/or  when  the 
denominator  for  the  spectrum  of  encounter  is  equal  to  zero  (spreading  function). 


Table  2.  Horizontal  Force  Comparison 


Horizontal  Force  Differential  (Connection-Ratio) 
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0.1 

-28.60% 

0.5 

28.60% 

1.0 
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0.1 

0.00% 

0.5 

20.00% 

1.0 

D) 

5 

0.1 
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0.5 

28.60% 

1.0 

'cD 

X 

10 

0.1 

-16.60% 

0.5 

16.60% 

1.0 

CD 

10 

0.1 

-9.10% 

0.5 

27.30% 

1.0 

cc 

> 

10 

0.1 

-14.30% 

0.5 

33.30% 

1.0 

> 

20 

0.1 

-20.00% 

0.5 

20.00% 

1.0 

CO 

o 

20 

0.1 

-12.00% 

0.5 

12.00% 

1.0 

H — 

'c 

20 

0.1 

-14.30% 

0.5 

33.30% 

1.0 

D) 

30 

0.1 

-12.50% 

0.5 

12.50% 

1.0 

30 

0.1 

-18.20% 

0.5 

18.20% 

1.0 

30 

0.1 

-14.30% 

0.5 

14.30% 

1.0 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  response  of  a  ship  advancing  in  a  stochastic  seaway  is  obviously  a 
complicated  phenomenon  involving  the  distinct  and  recognizable  interaction  between 
hydrodynamic  forces  and  vessel  dynamics.  Quantifying  and  measuring  either  sides  of  the 
equilibrium  equation  (39)  may  never  be  a  perfect  science,  however,  experiment  and 
prediction  prove  theoretical  insight  is  a  must  for  advances  in  naval  architecture.  While 
the  present  discussion  involves  only  distinct  evidence  for  vessel  response  to  waves  of 
various,  albeit  steady,  frequencies,  we  can  subsequently  still  predict  the  statistics  of  the 
responses  to  actual  random  seaways. 

The  measured  distribution  of  hydrodynamic  forces  along  the  length  of  the  vessels 
(strip  theory  predictions)  for  lower  frequencies  (longer  wavelengths)  become  greater  and 
more  adversely  affected  by  the  effects  of  three-dimensionality  and  forward  speed.  In  all 
cases,  lower  wave  periods  have  the  greatest  impact  on  the  vertical  and  horizontal 
connection  forces  (on  the  order  of  3.0).  For  higher  frequencies  the  agreement  between 
the  strip  theory  approximation  and  experimentation  are  reasonably  good  for  zero  forward 
speed,  but  the  agreement  between  theory  and  experiment  diverges  as  the  speed  increases 
(Beck). 

As  can  be  seen  by  the  speed  polar  plots,  forward  speed  dependence  is  clearly 
visible  in  the  numerical  results.  Pitch  motion  becomes  very  small  in  beam  seas,  but 
become  most  severe  in  head  seas.  The  effects  of  wave  frequency  and  forward  speed  in 
head  seas  have  been  illustrated.  In  general,  a  reduction  in  forward  speed  will  reduce  the 
heave  and  pitch  motion  in  long  wavelengths  and  increase  the  motion  in  short  waves.  The 
motions  in  short  waves  (L/L<1.0)  are  significantly  less  than  in  longer  wavelengths. 
Vertical  connection  forces  remained  relatively  the  same  for  increases  in  the  connection 
ratio  varying  only  a  maximum  of  9.1%.  Again,  keeping  all  parameters  constant  and 
varying  the  connection  ratio,  the  horizontal  components  on  the  other  hand  ranged  from  0- 
33%.  The  most  important  results/parameters  for  seakeeping  responses,  and  amplitude 
motbn  in  random  seas  are  wave  period,  wave  amplitude  (significant  wave  height),  and 
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forward  speed,  and  a  weaker  function  of  connection  ratio  for  the  ranges  considered. 
These  results  compare  with  results  found  earlier  (Okan),  whereby  the  connection  forces 
are  not  a  monotonic  function  of  the  towline  length  (therefore,  it  is  possible  to  select  a 
towline  length  so  that  the  connection  force  in  minimized). 

For  the  translational  modes,  the  surge  and  heave  motions  were  maximized  in  seas 
forward  of  the  bow.  Horizontal  plane  motions  tended  to  be  damped  out  while  vertical 
plane  motions  grew  larger  with  respect  to  speed  for  wave  angles  forward  of  beam.  For 
sway,  it  was  beam  seas  that  excited  the  largest  motions  while  yaw  was  greatest  in 
quartering  seas.  Sway  and  yaw  have  no  hydrostatic  restoring  forces  displaying  high 
amplitude  motions  in  quartering  and  beam  seas.  The  wave  encounter  frequencies  are  also 
much  lower  in  quartering  and  following  seas  than  in  head  seas  or  on  the  bow  and 
evidently  lead  to  larger  amplitude  sway  and  yaw  motions. 

It  must  be  remembered  that  the  incident  and  diffracted  wave  velocity  potential’s 
(39)  (respective  wave  velocity  potentials  before  and  after  contact  is  made  with  vessel) 
unsteady  components  are  all  higher  order  and  are  neglected  in  the  derivation  of  excitation 
forces  using  the  linear  theory  assumption.  Viscous  effects  also  play  a  more  important 
role  in  the  transverse  modes  than  in  the  horizontal.  For  example,  if  the  ship  is  at  an 
apparent  angle  of  attack  and  coincidently  yaws  relative  to  the  ship’s  forward  motion, 
viscous  effects  will  undoubtedly  come  into  effect  at  the  bow.  Significant  transverse 
rotation,  relative  to  angle  of  attach,  will  transpire.  However  such  motions  have  not  been 
incorporated  into  the  ship  motion  algorithm.  Additionally,  surge  motion  is  neglected  due 
to  the  simplified  head  seas  case  where  the  vessel  surges  forward  on  the  surface  of  the 
wave,  also  known  as  surfing.  Completing  the  free- body- diagram  on  the  vessel/wave 
interface,  it  is  easy  for  one  to  picture  the  buoyancy  force  normal  the  wave  face  and 
opposing  the  vessel  body  weight  component.  The  included  angle  between  the  normal 
force  and  the  weight  component  is  very  small,  none-the-less  the  included  positive  vector 
is  a  surge  force  and  is  clearly  present.  Sufficiently  large  tension  values,  however,  may 
stabilize  and  counter  such  forces  in  the  coupled  two -vessel  system 
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Continuing  with  our  discussion,  for  very  long  waves  the  heave  motion  amplitude 
approaches  the  wave  amplitude  and  the  pitch  amplitude  is  the  same  as  the  maximum 
wave  slope  for  such  conditions  where  the  vessel  moves  with  the  waves  (planar  motion) 
resulting  in  lower  transverse  forces.  As  the  ships  are  modeled  in-tow,  they  move  together 
similarly  at  the  corresponding  connection  points  as  a  result  of  the  coupling  connection. 
As  the  connection  ratio  decreases,  the  leading  and  trailing  ship  motions  become  closely 
matched  as  the  towline  length  is  reduced,  while  some  variation  is  observable  for  larger 
towline  lengths.  Furthermore,  the  horizontal  connection  force  becomes  larger  for  towing 
lengths  (1/L>0.1),  while  the  vertical  force  decreases  (of  varying  inconsistent  magnitude). 
Additionally,  there  is  no  clear  dependence  of  the  connection  force  and  consequent  ship 
motions  due  to  forward  speed,  there  is  a  large  variation  with  no  evident  trend. 

The  complexities  of  wave  resistance,  incident  waves  and  associated  singularities 
(wave  breaking),  as  well  as  some  of  the  nonlinearities  associated  with  the  free  surface 
problem  tend  to  make  the  speed  polar  plots  appear  more  as  a  glimpse  into  the  amplitude 
motions  of  six  degrees  rather  than  absolute.  While  the  most  advanced  present-day 
techniques  are  beyond  the  capacity  of  modern-day  computers,  and  will  likely  be  the  case 
for  years  to  come,  perhaps  the  simplex  method  describe  here  will  offer  a  numerical 
breach  into  the  many  unknowns  and  possibilities  of  ship  response  to  random  impulses.  In 
fact,  in  the  open  ocean,  the  evolution  and  development  of  waves  from  ripple  to  sea  state 
ten  is  likely  never  to  be  observed  given  the  obvious  explanation  that  the  wave  source  is, 
in  fact,  rarely  an  impulse.  Although  in  a  physical  sense  the  ocean  wave  and  its  multiple 
sea  states  cannot  be  reproduced  in  a  laboratory,  from  a  mathematical  standpoint  this  is 
quite  feasible.  “In  fact,  the  theoreticians  have  become  so  bold  as  a  result  of  the  success 
of  their  complicated  equations  that  there  is  danger  the  study  of  waves  will  fall  entirely 
into  the  hands  of  those  who  have  never  [even]  been  to  sea.”  (Bascom) 
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B.  RECOMMENDATIONS 


The  overall  effect  of  vessel  control  must  be  realized  in  order  to  better  determine 
the  surge  and  roll  motion  amplitudes.  Surge  is  severely  affected  by  propeller  thrust, 
rudder  design,  and  steering  gear  response.  Propeller  response  and  the  unsteady  wake 
behind  SLICE  would  either  serve  to  dampen  or  excite  mode  amplitudes  for  the 
KAIMALINO.  Furthermore,  the  idealized  mechanical  linkage  between  the  two  vessels 
is  a  major  contributor  to  the  coupling  mechanics  between  the  vessels  and  correlates 
directly  to  roll  motion.  The  damping  arising  from  wave  making  due  to  roll  is  very  small 
for  this  mode;  while  for  heave,  yaw,  sway,  etc.  it  is  a  major  contributor.  This  idealized 
linkage  assumption,  however,  should  be  altered  to  incorporate  a  more  realistic 
damping/excitation  force  when  considering  the  roll  motion. 

Numerically,  the  boundary  conditions  on  the  hull  and  on  the  free  surface  both 
contained  terms  based  on  the  steady  perturbation  potential.  In  trying  to  solve  this 
boundary  value  problem,  it  is  determined  that  there  exists  an  interaction  between  the 
steady  and  unsteady  components  of  this  potential.  Addressing  the  nonlinearities 
associated  with  the  potential  function,  in  particular,  the  velocity  squared  term  in  the 
pressure/Bernoulli  equation  and  viscous  (damping)  effects.  Similarly,  the  smooth  free 
surface  assumption  should  be  reconsidered  (incident  wave  brake).  Unique  solutions  of 
the  Laplace  equation  (Appendix  B)  require  boundary  conditions  on  all  surfaces 
surrounding  the  fluid  domain.  The  major  nonlinearities  in  the  general  potential  flow 
problem  are  in  the  free- surface  boundary  conditions  that  involve  the  square  of  the  fluid 
velocities  and  products  of  these  velocity  terms  with  the  unknown  free- surface  amplitude. 
This  potential  flow  problem  still  is  quite  difficult  and  further  underlying  assumptions  will 
be  necessary.  For  our  simplifying  assumption  of  linear  theory,  the  interaction  was 
considered  negligible  and  only  the  free-stream  value  was  deemed  necessary.  Although 
many  in  the  hydrodynamics  field  practice  this  assumption;  is  this  a  reasonable  hypothesis 
considering  our  significant  wave  heights  and  varying  ship  velocities? 

Any  disturbance  of  the  water  surface,  including  the  passage  of  a  ship,  produces 
waves.  Much  of  the  power  expended  in  propelling  the  ship  goes  into  such  said  wave- 
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making  processes  and  while  this  is  a  prime  concern  of  the  marine  engineer  in  the  interest 
of  inefficiencies  and  power  losses,  it  may  also  be  of  interest  to  future  proceedings  and 
research  into  using  SWATH  vessels  in- tow.  As  the  ship  passes  through  the  water  it  is 
accompanied  by  at  least  three  different  pressure  disturbances  on  each  side,  which  produce 
several  distinct  wave  trains  in  consistent  order,  direction,  and  symmetry.  Lord  Kelvin 
investigated  the  pattern  of  waves  generated  by  an  idealized  pressure  disturbance 
concentrated  at  a  point  and  moving  in  a  straight  line.  The  Kelvin  wave  pattern,  as  it  is 
called,  is  characterized  by  (1)  diverging  waves  (a  series  of  curved  crests,  concave 
outward  and  lying  in  echelon  position);  (2)  transverse  waves  (convex  forward  and 
perpendicular  to  the  direction  of  motion);  (3)  crest  intersections  (where  the  diverging  and 
transverse  waves  coincide).  It  is  this  trailing  pattern  that  should  be  further  researched  and 
incorporated  into  the  two- vessel  design.  The  Volex  wave  system  (ship  applications  as 
apposed  to  point  disturbance)  and  pressure  deductions  should  be  determined  for  the 
SLICE  twin  hull  configuration.  Do  the  resultant  stern  waves  contribute  to  the  overall 
wave  resistance?  When  considering  the  positive  pressure  wave  at  the  stem,  should  this 
attribute  to  added  resistance  and  be  incorporated  into  the  connection  force? 

Finally,  should  a  more  realistic  six  degree- of- freedom  coupling  be  developed,  the 
potential  to  model  this  result  at  the  Fbet  Numerical  Oceanographic  Center  (FNOC), 
Monterey,  CA,  in  real-time  has  been  addressed  to  FNOC,  and  may  better  ascertain  the 
feasibility  and  validity  of  the  coupling  algorithm.  (Virtually  place  the 
SFICE/KAIMAFINO  coupling  anywhere  in  the  world.) 
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Anecdotal: 


With  the  advantages  of  the  computational  domain  come  the  disadvantages.  It  has 
been  found,  without  careful  observation  numerical  waves  can  reflect  from  the  boundary 
conditions  and  return  to  and  interfere  with  initial  conditions  fouling  results  (Beck  and 
Reed,  2001). 

“To  impact  ship  design,  computational  methods  must  be  integrated  into  the  design 
process.  This  in  turn  means  that  the  computer  codes  must  be  fast  enough  that  they  can 
provide  results  in  no  more  than  an  hour  or  two  of  computing  time.  While  there  has  been 
work  done  on  improving  the  computational  efficiency  of  computer  codes  for  ship 
hydrodynamics  problems,  a  tremendous  amount  of  work  still  needs  to  be  done.  The  work 
to  date  on  computational  efficiency  has  resulted  in  a  large  enough  literature  base  that  this 
subject  is  worthy  of  a  paper  of  its  own;  time  and  space  [and  money]  obviously  preclude 
this.”  (Beck  and  Reed,  2001). 
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APPENDIX  A.  SPEED  POLAR  PLOT  RESULTS 


Considering  the  limited  memory  presently  alloeated  to  MS  WORD,  only  a  small 
sampling  (20  modal  solutions)  of  the  360  member  set  will  be  ineluded  to  illustrate  the 
numerieal  results.  If  further  detail  is  required  of  the  results  outlined  in  Table  1,  please 
email  Professor  Fotis  Papoulias  (papoulias@nps.navv.mil),  Meehanieal  Engineering 
Department,  Naval  Postgraduate  Sehool;  Monterey,  CA. 
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Figure  A- 1.  Vertical  Connection  Force  (HS:  5  feet,  Tn,:  5  sec.,  1/L:  0.1) 


Figure  A-2.  Transverse  Connection  Force  (HS:  5  feet,  Tmi  5  sec.,  1/L:  0.1) 
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Leading  Ship,  Heave 


Figure  A-3.  SLICE  Heave  Amplitude  Motion  (HS:  5  feet,  T^:  5  sec.,  1/L:  0.1) 
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Figure  A-4.  KAIMALINO  Heave  Amplitude  Motion  (HS:  5  feet,  Tn,:  5  sec.,  1/L: 

0.1) 
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Leading  Ship,  Pitch 


Figure  A-5.  SLICE  Pitch  Amplitude  Motion  (HS:  5  feet,  Tn,:  5  sec.,  1/L:  0 


Trailing  Ship,  Pitch 


Figure  A-6.  KATMATJNO  Pitch  Amplitude  Motion  (HS:  5  feet,  Tn,:  5  sec., 
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Leading  Ship,  Sway 
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Figure  A-7.  SLICE  Sway  Amplitude  Motion  (HS:  5  feet,  Tn,:  5  sec.,  1/L:  0.1) 
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Figure  A-8.  KAIMALINO  Sway  Amplitude  Motion  (HS:  5  feet,  Tmi  5  sec.,  1/L: 

0.1) 
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Leading  Ship,  Yaw 


Figure  A-9.  SLICE  Yaw  Amplitude  Motion  (HS:  5  feet,  Tn,:  5  sec.,  1/L:  0.1) 

Trailing  Ship,  Yaw 


Figure  A- 10.  KAIMALINO  Yaw  Amplitude  Motion  (HS:  5  feet,  Tmi  5  sec.,  1/L: 

0.1) 
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Vertical  Connection  Force 
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Figure  A-351.  Vertical  Connection  Force  (HS:  30  feet,  Tm)  20  sec.,  1/L:  1.0) 
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Figure  A-352.  Transverse  Connection  Force  (HS:  30  feet,  Tn,:  20  sec.,  1/L:  1.0) 
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Figure  A-353.  SLICE  Heave  Amplitude  Motion  (HS:  30  feet,  T^:  20  sec.,  1/L:  1.0) 
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Figure  A-354.  KAIMALINO  Heave  Amplitude  Motion  (HS:  30  feet,  Tm:  20  sec.,  1/L: 

1.0) 


36 


Figure  A-355.  SLICE  Pitch  Amplitude  Motion  (HS:  30  feet,  Tn,:  20  sec.,  1/L:  1.0) 
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Figure  A- 356.  KAIMALINO  Pitch  Amplitude  Motion  (HS:  30  feet,  Tn,:  20  sec.,  1/L: 

1.0) 
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Leading  Ship,  Sway 


Figure  A-357.  SLICE  Sway  Amplitude  Motion  (HS:  30  feet,  T„,:  20  sec.,  1/L:  1.0) 


Trailing  Ship,  Sway 


Figure  A-358.  KAIMALINO  Sway  Amplitude  Motion  (HS:  30  feet,  Tn,:  20  sec.,  1/L: 

1.0) 
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Leading  Ship,  Yaw 


Figure  A-359.  SLICE  Yaw  Amplitude  Motion  (HS:  30  feet,  Tn,:  20  sec.,  1/L:  1.0) 
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Figure  A- 360.  KAIMALINO  Yaw  Amplitude  Motion  (HS:  30  feet,  Tmi  20  sec.,  1/L: 

1.0) 
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APPENDIX  B.  FUNDAMENTAL  HYDRODYNAMIC  THEORY 


The  six-degree  of  freedom  motion  of  the  SLICE/KAIMALINO  vessel 
configuration  in- tow  requires  a  fairly  good  understanding  of  the  three  dimensional 
motion  and  the  corresponding  wave  hydrodynamics.  In  order  to  investigate  the  motions 
on  the  ocean  surface  (particularly  wave  velocity),  we  must  first  develop  the  basic 
equations  (and  understanding)  of  the  hydrodynamics  of  an  ideal  (nonviscous)  fluid.  In 
our  study,  we  assumed  that  the  water  is  incompressible,  the  density  is  constant  and  the 
flow  is  irrotational.  These  are  standard  assumptions  in  seakeeping  studies.  We  also 
assumed  that  the  ships  operate  in  infinitely  deep  water.  The  characterizations  are  as 
follows:  (1)  kinematic  equations;  (2)  conservation  of  mass  (continuity  equation);  (3) 
dynamic  equations  or  Bernoulli’s  equation  describing  irrotational  motion  of  an 
incompressible  fluid. 

Using  lagrangian  techniques,  we  can  describe  any  number  of  particles  and  their 
relative  positions  with  respect  to  the  fluid  motion.  Letting  E  =  e(v,  y,z,t)  stand  for  a 
(scalar)  property  of  the  ocean  wave  medium,  the  lagrangian  rate  of  change  for  a  particle 
function,  E,  for  some  fixed  point  ro: 

=  EU„,y„,z„,tJ  (1) 


When  a  particle  at  point  r  at  the  time,  t,  arrives  at  a  neighboring  point  at  the  time, 
t-i-dt,  the  function  E  is  the  following: 


^[x^  +  +  +  ,t^  +  dt)  =  E^  + 


dl  + 


Jo 


dt]/  -I- 


\dz  / 


dC  + 


dt  (2) 


Whereby  the  time  rate  of  change  of  E  from  the  lagrangian  point  of  view: 


dE  dE 

—  =v  •  VE-h — 
dt  dt 


(3) 


where. 
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(4) 


.dx  .  3vi/  , 

at  at  at 

.aE  .aE  ,  aE 

ax  ay  az 


Thus  is  described  the  general  kinematic  equation  of  a  water  particle. 

For  our  assumptions  of  an  incompressible  fluid,  which  is  an  accurate  model  of  the 
ocean  medium  at  the  surface  where  the  pressure  changes  are  minimal,  the  continuity 
equation  resolves  to: 

V-v=0  (5) 

where  v  is  the  displacement  velocity  which  is  derivable  from  the  velocity 
potential,  ([)  ,  whereby  the  velocity  potential  in  an  incompressible  fluid  must  satisfy 
Laplace’s  equation: 

V  ■  V4)  =  0  (6) 

The  surface  wave  is  the  manifestation  of  pressure  changes  and  water- particle 
motions  affecting  the  entire  body  of  fluid-  in  a  column  sense  all  the  way  through  its 
respective  depth,  surface  to  ocean  floor.  The  motion  of  particles  under  the  idealized 
given  conditions  can  be  characterized  by  this  velocity  potential,  (|)  .  From  this  function 
all  desired  wave  properties  can  be  determined. 

We  shall  now  derive  the  Bernoulli  Equation  by  applying  Newton’s  second  law  to 
the  mass  pAx  in  an  element  of  volume  Ax  of  a  fluid.  The  net  force  acting  on  the  element 
to  give  it  an  acceleration  dv/dt  includes  an  external  body  force,  Fe  (due  to  gravity),  and  an 
internal  force  Fj  (pressure  gradients  within  the  fluid).  The  pressure  force  between  each 
element  can  be  described  as  follows: 


AF.  =  -VP  Ax  (7) 

AF^  =  -  pVQAx  (8) 

where  Q.  is  the  gravitational  potential  energy  per  unit  mass,  g=-Vf2  (i.e.  positive 
y-axis  is  directed  upward,  then  f2=gy,  and  g=-jg. 
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F=ma  (9) 


-  yPAx  -  pVf2Ax  =  p  —Ax 

dt 


(10) 


referring  to  our  earlier  discussion  of  equation  (3), 


dE  _  dE 
—  =  V  •  VE  +  — 
dt  dt 


(11) 


we  may  now  express  the  acceleration  as  follows: 

dv 

-VP -  pVQ.- V  Vv  +—  (12) 


When  generally  discussing  waves  on  the  surface  of  a  liquid,  it  is  customary  to 
assume  that  the  liquid  particle  motion  is  irrotational.  This  means  that  the  curl  of  V 
everywhere  vanishes,  and  essentially  the  mathematical  characteristics  of  eddies,  or 
vorticites  is  lost  and  obviously  not  of  our  concern.  When  described  in  an  ideal  fluid, 
these  vortices  persist  forever  to  conserve  angular  momentum;  however,  in  reality  vortices 
die  out  due  to  fluid  viscosity  (i.e.  smoke  rings).  Therefore, 

Vxv=0  (13) 

and, 

(Vxv)xv  =v -Vv  -  Vv -v  =0  (14) 
so  that: 

V  -  Vv  =  (Vv)-v  =  V(|v^)  (15) 

The  no  curl  condition  is  a  necessary  and  sufficient  condition  that  v  can  be  derived  from 
the  velocity  potential  (|). 


V  =  -V^  (16) 

and  since. 


dt 


dt 


(17) 


equation  10  becomes: 
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VP-Vi^  =  -V^+v(Yv")  (18) 
p  at 

integration  can  then  be  carried  through  by  dotting  dr  into  it  and  noting  that: 
dE  =  dr-VE  (19) 

where  dE  is  the  difference  between  two  neighboring  points  separated  by  dr . 
Equation  18  becomes: 

-—dP-da  =  -d^  +  d{iv^)  (19) 
p  at 

and  when  integrated  gives  the  relation: 

(dP  d(|)  ,  , 

=  (20) 

p  at 

and  when  considered  to  be  incompressible,  we  resolve  to  Bernoulli’s  Equation: 

P  1  , 

-+a+W  =  ^+F(t)  (21) 

p  dt 

whereby  the  constant  of  integration,  and  any  time- dependent  component  has  been 

3(|) 

absorbed  into  the  undetermined  —  component.  Taking  the  time  derivative  and  treating 

at 

P  as  a  constant  and  considering: 


V. 


dt  dx 

3\|/  34) 


V  =  —  =  - —  (22) 
dt  dy  ^  ’ 


V,  = 


3C  3(^ 


dt  dz 

evaluating  y=h,  we  find: 

34)  1  3^4^ 


(23) 


dy  g  dr 

showing  that  velocity  depends  on  position  and  time,  Bernoulli’s  equation  can  now 


be  used  to  compute  pressure  at  any  point  in  the  fluid.  Solving  for  the  velocity  potential 
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function  and  knowing  that  it  is  a  function  of  position  and  time,  we  must  satisfy  the 
Laplace  equation: 


3^(1)  3^(|) 

dx^  ^  dy^ 


=  0  (24) 


where; 

(^Uy,t)=X(^)-F(0-r(0  (25) 


The  spatial  functions  are  determined  from  separation  of  variables,  and  the  general 
solution: 


X{x)  = 

F(y)  =  Ce'^+De-'^ 


(26) 


considering  the  boundary  conditions  at  y=0: 
3(b  dY 

which  requires  that  C=D,  where; 

F(y)=  2Ccosh(Ky)  (28) 


similarly  the  velocity  potential  has  the  form: 

4)  =  Acosh  Kye'“r(0  (29) 

and  substituting  into  equation  24,  we  obtain  the  equation  for  T(t): 
d  ^7^ 

— 5-+(gKtanh  idz)r  =  0  (30) 
dt 

which  leads  to  the  simple  solution  of  harmonic  time  dependence  with  the  angular 
frequency: 

a)  =  (gKtanhidz)'/"  (31) 
wave  velocity  is  then: 
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— tanh  Kh 


]/2 


(32) 


for  the  case  of  deep  water,  say  when  h>X/2,  tanhKh>tanh7C  =0.996,  so  that  the 
expression  for  the  wave  velocity  becomes  rather  closely: 


,(X<h)  (33) 
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APPENDIX  C.  VERTICAL  PLANE  COUPLING  DERIVATION 


Defining  the  equations  of  motions  for  a  point  on  the  each  of  the  SLICE  and  the 
KAIMALINO: 

[A,+A,]{fi-}  +  [5,]{fi}  +  [CJ{ri}  =  {E,J+{/,}  (1) 

[A,  +  A^]{fi}  +  [B^m  +  [C,]{ri}  =  {FJ  +  {f^}  (2) 

where  /^and  are  the  connection  forces  acting  on  SLICE  and  KAIMALINO 
respectively,  where  =  -/^ .  The  connection  forces  and  absolute  motions  are 
illustrated  in  Eigure  C  1 . 


Figure  C  1.  Connection  Force  and  Motions  of  the  Ships  (From:  Orhan) 


Continuing  with  the  derivation  of  Equation  (39): 

[^^\A+A)+m^B  +  Cm  =  {FJ  +  [f}  (3) 

Since  the  exponential  exists  in  all  terms,  it  is  canceled  and  the  equations  of 
motion  in  the  frequency  domain  become  as  follows: 

A  =  [-(of(A  +  A)  +  m^B+C]  (4) 
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We  know  that  A  is  a  6x6  matrix  for  six  degrees  of  freedom  system: 

Ati  =  {F}  (5) 

tl  =  mv(A){F}  (6) 

where; 

{F}  =  {FJ  +  {f}  (7) 

Ship  motions  due  to  regular  waves  of  a  given  wavelength  and  direction  are  now 
determined  for  a  given  forward  speed  (V).  Motions  in  vertical  and  horizontal  planes  are 
decoupled  for  our  discussion  and  may  be  solved  as  two  distinct  3x3  systems.  The  heave 
and  pitch  modes  and  the  resultant  interactions  will  be  derived  from  this  point  forward 
(sway  and  yaw  are  similar).  The  expanded  equations  of  motion  in  two  degrees  of 
freedom  become  as: 


A33 

As 

4' 

■  A+/' 

_Ag3 

4. 

As. 

fs  +  A. 

Now  we  can  define  the  equations  of  motion  in  vertical  plane  (two  degrees  of 
freedom)  for  the  two  vessels  as: 

A33,.ri3^5  +A35,sri5^  =  7^3  ^  +  fs  (12) 

A53,sTl3  ^  +A55,5Tl5  ^  —  T’j  ^  ~  f S  (^3) 

A33.xri  3  ^  +  A35,xri  5  ^  =  7^3  ^  +  /^  (14) 

AsiatT]  3  ^  +  Ass.xTIj  3,  —  F^jr  —  (15) 
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where  x  is  the  distance  between  the  connection  point  and  the  center  of  gravity  of 
that  vehicle.  In  order  to  solve  these  four  equations,  we  make  the  following  substitution: 


^3,S  “ 

M3,J 

+v. 

3,5 /s 

^5,S  “ 

M5,S 

+v. 

5,5 /j 

^3, A"  ~ 

M3,a: 

+v 

3,a:/a' 

= 

M5,a: 

+v 

5,a:/a' 

where; 

(16) 


fs=-fK=-f  (17) 

Using  Cramer’s  Rule,  the  absolute  motions  can  be  described  by  (Nash): 
^5  =Tl3,S 

-  -  (lo) 

“^3, a: 

When  we  combine  (13)  and  (14)  as: 

~  ^3,S  3,S  f  ~  ^5,S^S  5,S^S  f 

~  M'S, a:  3,A:/  ~  \^5,K^K  5,K^Kf 

Let  a  in  terms  of  p,  ’s  and  b  in  terms  of  V ’s  be  defined  as: 

^  ~  Ms.s  ~  Ms, a:  M5,A'^A'  ~  M5,S-^S 

^  ~  3,S  3, K  5,K  5,S 


(20) 


Therefore,  our  equation  of  motion  becomes: 

^s-^K=ci-bf  (21) 

Again  referring  to  Figure  2, 

/  ^T^s-^K  (22) 
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Combining  (17)  and  (18)  gives  us  our  final  equation,  which  gives  us  the 
connection  force  between  two  ships: 


/=7^  (23) 

-+b 

T 


Figure  C  2.  KAIMALINO  Tension  vs.  Speed  Plot  (After:  Lockheed  Martin) 
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APPENDIX  D.  MATLAB  COUPLING  ALGORITHM 


%  Horizontal  Plane 
%  Dimensional  version  (U.S.  units) 

%  Contour  plots  (heading/speed) 

%  Two  parameter  Bretschneider  spectrum  -  Short  crested  seas 


% 

%  Get  run  info 
% 

HS=input('Significant  Wave  Height  (feet)  =  '); 

T_m=input('Modal  Period 
l=input('Length  (l/L)  =  '); 

% 

omega_m=2*pi/T_m; 

HS_string=num2str(HS); 
warning  off 
% 

lambda_min=20; 
lambda_max=1000; 
delta_lambda=20; 
rho=1.9905; 
zeta=l; 

L=105; 
g=32.2; 
x_s=-46; 
x_k=+40; 
beta_incr=5; 
motion_ratio=2; 

% 

%  The  matdata  output  files  default  to  the  vertical  only  format  when  the 
%  heading  angle  is  0  or  180  degrees. 

%  Set  up  file  reading  format. 

% 

% 

%  GENERAL  DATA 

% 

lambda=lambda_min:delta_lambda:lambda_max;  %  Vector  of  wavelengths 

wavenumber=2.0*pi. /lambda;  %  Wave  number 

omega=sqrt(wavenumber*g);  %Wave  frequency 

period=2.0*pi. /omega; 

omega=omega'; 

filesize=size(lambda); 

if  beta_incr=  =  15 

load  data_15;  %  Load  main  data  file,  15  deg.  heading  increments 

beta_incr_values=13; 
end 

if  beta_incr==5 

load  data_5;  %  Load  main  data  file,  5  deg.  heading  increments 

beta_incr_values=37; 
end 
% 


%  Min  wave  length  (ft) 

%  Max  wave  length  (ft) 

%  Wave  length  increment  (ft) 

%  Water  density 
%  Regular  wave  height 
%  Reference  length  for  nondimensionalization 
%  Gravitational  constant 
%  FRONT  SHIP  attachment  point 
%  REAR  SHIP  attachment  point 
%  Increment  in  sea  direction  (deg) 

%  Ratio  for  motion  plotting  comparison 
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%  Random  wave  calculations 

%  Bretschneider  spectrum  -  weighted  by  the  spreading  function 
% 

iSpeed=0; 

for  V_knots=0:l:20,  %  Loop  on  speed 

iSpeed  =  iSpeed  +  l; 

% 

%  Get  tension  from  curvefitting  data. 

%  Applicable  for  speeds  between  1  and  20  ft/sec. 

% 

T=-1.762*V_knots"4+63.675*V_knots"3-580.8*V_knots"2  +  2485.9*V_knots-34.047; 
% 

V_string=num2str(V_knots) 

V=V_knots*  1.6878;  %  Convert  to  ft/sec 

A=(1.25/4)*(omega_m^4)*(HS"2); 

B=1.25*omega_m^4; 

S_main=(A./omega.^5).*exp(-B./omega.^4); 

% 

ibeta=0; 

for  beta=0:beta_incr:360,  %  Loop  on  sea  direction 

beta_set=beta; 
if  beta>180 

beta_set=360-beta; 

end 

ibeta=ibeta+l; 

beta; 

% 

%  Spreading  between  (beta-90)  and  (beta-l-90) 

% 

beta_spread_low  =beta_set  -  90; 
beta_spread_high=beta_set  -l-  90; 
ibeta_spread=0; 

% 

%  Loop  between  (beta-90)  and  (beta-l-90) 

%  Use  either  15  deg  increments  (13  values  total)  or 
%  5  deg  increments  (total  of  37  values) 

% 

for  beta_spread=beta_spread_low:beta_incr:beta_spread_high, 
ibeta_spread=ibeta_spread  -i-  1; 
beta_spread_read=beta_spread; 
beta_spread_vector(ibeta_spread)=beta_spread; 
if  beta_spread<0 

beta_spread_read=-beta_spread; 

end 

if  beta_spread>180 

beta_spread_read=360-beta_spread; 

end 

beta_string=num2str(beta_spread_read); 

trigg=30; 

f2loc=26;  f6loc=30; 
if  beta_spread_read==0 
trigg=27; 

f2loc=25;  f6loc=27; 
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elseif  beta_spread_read==180 
trigg=27; 

f2loc=25;  f6loc=27; 
end 

lambda_size=trigg*filesize(2); 

% 

%  Load  FRONT  SHIP  data  file  msvhV_beta.txt 
% 

load_filename=strcat('msvh',V_string,'_',beta_string); 

filename_s=eval(load_filename); 

% 

%  Load  REAR  SHIP  data  file 

% 

load_filename=strcat('mkvh',V_string,'_',beta_string); 

filename_k=eval(load_filename); 

omegae=omega'-wavenumber*V*cos(beta_spread*pi/180);  %Freq.  of  Encounter 

omegae=omegae'; 

periode=2.0*pi./omegae; 

% 

%  HORIZONTAL  PLANE  RESPONSE  CALCULATIONS 

% 

%  SLICE 

% 

%  Set  mass  matrix  elements 
% 

M22s=filename_s(2:trigg:lambda_size,2); 

M26s=filename_s(2:trigg:lambda_size,6); 

M62s=filename_s(6:trigg:lambda_size,2); 

M66s=filename_s(6:trigg:lambda_size,6); 

% 

%  Added  mass  terms 
% 

A22s=filename_s(8:trigg:lambda_size,2); 

A26s=filename_s(8:trigg:lambda_size,6); 

A62s=filename_s(12:trigg:lambda_size,2); 

A66s=filename_s(12:trigg:lambda_size,6); 

% 

%  Damping  terms 
% 

B22s=filename_s(14:trigg:lambda_size,2); 

B26s=filename_s(14:trigg:lambda_size,6); 

B62s=filename_s(18:trigg:lambda_size,2); 

B66s=filename_s(18:trigg:lambda_size,6); 

% 

%  Hydrostatic  terms 
% 

C22s=filename_s(20:trigg:lambda_size,2); 

C26s=filename_s(20:trigg:lambda_size,6); 

C62s=filename_s(24:trigg:lambda_size,2); 

C66s=filename_s(24:trigg:lambda_size,6); 

% 

if  beta_spread==0 

F2s_t=zeros(50,l);  F6s_t=zeros(50,l); 
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elseif  beta_spread==180 

F2s_t=zeros(50,l);  F6s_t=zeros(50,l); 
else 
% 

%  Total  exciting  forces 
% 

F2s_t_amp=fllename_s(f2loc:trlgg:lambcla_slze,5); 

F6s_t_amp=fllename_s(f6loc:trlgg:lambcla_slze,5); 

F2s_t_pha=fllename_s(f2loc:trlgg:lambcla_slze,6); 

F6s_t_pha=fllename_s(f6loc:trlgg:lambcla_slze,6); 

F2s_t=F2s_t_amp.*exp(i*F2s_t_pha.*pl/180.0); 

F6s_t=F6s_t_amp.*exp(i*F6s_t_pha.*pl/180.0); 

% 

%  Froude/Krylov  exciting  forces 
% 

F2s_f_amp=fllename_s(f2loc:trlgg:lambda_slze,l); 

F6s_f_amp=fllename_s(f6loc:trlgg:lambda_slze,l); 

F2s_f_pha=fllename_s(f2loc:trlgg:lambda_slze,2); 

F6s_f_pha=fllename_s(f6loc:trlgg:lambda_slze,2); 

F2s_f=F2s_f_amp.*exp(i*F2s_f_pha.*pl/180.0); 

F6s_f=F6s_f_amp.*exp(i*F6s_f_pha.*pl/180.0); 

% 

%  Diffraction  exciting  forces 
% 

F2s_d_amp=filename_s(f2loc:trigg:lambda_size,3); 

F6s_d_amp=filename_s(f6loc:trigg:lambda_size,3); 

F2s_d_pha=filename_s(f2loc:trigg:lambda_size,4); 

F6s_d_pha=filename_s(f6loc:trigg:lambda_size,4); 

F2s_d=F2s_d_amp.*exp(i*F2s_d_pha.*pi/180.0); 

F6s_d=F6s_d_amp.*exp(i*F6s_d_pha.*pi/180.0); 

% 

end 

%  KAIMALINO 

% 

%  Set  mass  matrix  elements 
% 

M22k=filename_k(2:trigg:lambda_size,2); 

M26k=filename_k(2:trigg:lambda_size,6); 

M62k=filename_k(6:trigg:lambda_size,2); 

M66k=filename_k(6:trigg:lambda_size,6); 

% 

%  Added  mass  terms 
% 

A22k=filename_k(8:trigg:lambda_size,2); 

A26k=filename_k(8:trigg:lambda_size,6); 

A62k=filename_k(12:trigg:lambda_size,2); 

A66k=filename_k(12:trigg:lambda_size,6); 

% 

%  Damping  terms 
% 

B22k=filename_k(14:trigg:lambda_size,2); 

B26k=filename_k(14:trigg:lambda_size,6); 

B62k=filename_k(18:trigg:lambda_size,2); 
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B66k=filename_k(18:trigg:lambda_size,6); 

% 

%  Hydrostatic  terms 
% 

C22k=filename_k(20:trigg:lambda_size,2); 
C26k=filename_k(20:trigg:lambda_size,6); 
C62k=filename_k(24:trigg:lambda_size,2); 
C66k=filename_k(24:trigg:lambda_size,6); 
if  beta_spread==0 

F2k_t=zeros(50,l);  F6k_t=zeros(50,l); 
eiseif  beta_spread==180 

F2k_t=zeros(50,l);  F6k_t=zeros(50,l); 
eise 
% 

%  Totai  exciting  forces 
% 

F2k_t_amp=fiiename_k(f2ioc:trigg:iambda_size,5); 
F6k_t_amp=fiiename_k(f6loc:trigg:iambda_size,5); 
F2k_t_pha=fiiename_k(f2ioc:trigg:iambda_size,6); 
F6k_t_pha=fiiename_k(f6ioc:trigg:iambda_size,6); 
F2k_t=F2k_t_amp.*exp(i*F2k_t_pha.*  pi/ 180.0); 
F6k_t=F6k_t_amp.*exp(i*F6k_t_pha.*  pi/ 180.0); 

% 

%  Froude/Krylov  exciting  forces 
% 

F2k_f_amp=fiiename_k(f2loc:trigg:iambda_size,l); 
F6k_f_amp=fiiename_k(f6ioc:trigg:lambda_size,l); 
F2k_f_pha=fiiename_k(f2loc:trigg:iambda_size,2); 
F6k_f_pha=fiiename_k(f6ioc:trigg:lambda_size,2); 
F2k_f=F2k_f_amp.*exp(i*F2k_f_pha.*  pi/ 180.0); 
F6k_f=F6k_f_amp.*exp(i*F6k_f_pha.*  pi/ 180.0); 

% 

%  Diffraction  exciting  forces 
% 

F2k_d_amp=fiiename_k(f2ioc:trigg:iambda_size,3); 
F6k_d_amp=fiiename_k(f6ioc:trigg:lambda_size,3); 
F2k_d_pha=fiiename_k(f2loc:trigg:lambda_size,4); 
F6k_d_pha=fiiename_k(f6ioc:trigg:lambda_size,4); 
F2k_d=F2k_d_amp.*exp(i*F2k_d_pha.*  pi/ 180.0); 
F6k_d=F6k_d_amp.*exp(i*F6k_d_pha.*  pi/ 180.0); 
end 
% 

%  MATCHING  CONDITION 

% 

A22bar_s=-(omegae.^2).*(M22s+A22s)+i*omegae.*B22s+C22s; 

A26bar_s=-(omegae.^2).*(M26s+A26s)+i*omegae.*B26s+C26s; 

A62bar_s=-(omegae.^2).*(M62s+A62s)+i*omegae.*B62s+C62s; 

A66bar_s=-(omegae.^2).*(M66s+A66s)+i*omegae.*B66s+C66s; 

A22bar_k=-(omegae.^2).*(M22k+A22k)+i*omegae.*B22k4G22k 

A26bar_k=-(omegae.^2).*(M26k+A26k)+i*omegae.*B26k+C26k 

A62bar_k=-(omegae.^2).*(M62k+A62k)+i*omegae.*B62k+C62k 

A66bar_k=-(omegae.^2).*(M66k+A66k)+i*omegae.*B66k+C66k 

% 
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A26_s_den=(A22bar_s.*A66bar_s-A26bar_s.*A62bar_s); 

A26_k_den  =  (A22bar_k.*A66bar_k-A62bar_k.*A26bar_k); 

% 

mu2_s=(A66bar_s.*F2s_t-A26bar_s.*F6s_t)./A26_s_den; 

nu2_s=(A66bar_s-A26bar_s*x_s)./A26_s_den; 

mu6_s=(A22bar_s.*F6s_t-A62bar_s.*F2s_t)./A26_s_den; 

nu6_s=(A22bar_s*x_s-A62bar_s)./A26_s_den; 

mu2_k=(A66bar_k.*F2k_t-A26bar_k.*F6k_t)./A26_k_den; 

nu2_k=(A66bar_k-A26bar_k*x_k)./A26_k_den; 

mu6_k=(A22bar_k.*F6k_t-A62bar_k.*F2k_t)./A26_k_den; 

nu6_k=(A22bar_k*x_k-A62bar_k)./A26_k_den; 

% 

a=mu2_s+mu6_s*x_s-mu2_k-mu6_k*x_k; 

b=nu2_s+nu6_s*x_s+nu2_k+nu6_k*x_k; 

f=a./(l/T+b); 

% 

f_s=-f;  %  Connection  force  on  SLICE 

f_k=f;  %  Connection  force  on  KAIMALINO 

eta2_s=mu2_s+nu2_s.*f_s;  %  SLICE  sway 

eta6_s=mu6_s+nu6_s.*f_s;  %  SLICE  yaw 

eta2_k=mu2_k+nu2_k.*f_k;  %  KAIMALINO  sway 

eta6_k=mu6_k+nu6_k.*f_k;  %  KAIMALINO  yaw 

xi_s=eta2_s+eta6_s*x_s;  %  SLICE  motion  at  connection 

xi_k=eta2_k+eta6_k*x_k;  %  KAI  MALI  NO  motion  at  connection 

xi0_s=mu2_s+mu6_s*x_s;  %  SLICE  motion  at  connection  for  zero  f 

xi0_k=mu6_k+mu6_k*x_k;  %  KAIMALINO  motion  at  connection  for  zero  f 

deita_beta(ibeta_spread)=beta_set-beta_spread; 

S=S_main  *(2/pi)*(cos((beta_set-beta_spread)*pi/180))^2;  %  Muit.by  (2/pi)cos^2 

Se=S./abs((l-(2.0/g)*omega*V*cos(beta_spread*pi/180)));  %  S(w)  to  S(we) 

% 

%  Define  response  spectra  and  store  as  functions  of  theta 
% 

Sf(:,ibeta_spread)=((abs(f)).^2).*Se; 

Sxi_s(:,ibeta_spread)  =  ((abs(xi_s)).^2).*Se; 
Sxi_k(:,ibeta_spread)=((abs(xi_k)).^2).*Se; 
SxiO_s(:,ibeta_spread)=((abs(xiO_s)).^2).*Se; 

SxiO_k(:,ibeta_spread)  =  ((abs(xiO_k)).^2).*Se; 

SF2s_t(:,ibeta_spread)  =  ((abs(F2s_t)).^2).*Se; 

SF2k_t(:,ibeta_spread)  =  ((abs(F2k_t)).^2).*Se; 

Seta2_s(:,ibeta_spread)  =  ((abs(eta2_s)).^2).*Se; 
Smu2_s(:,ibeta_spread)=((abs(mu2_s)).^2).*Se; 

Seta6_s(:,ibeta_spread)  =  ((abs(eta6_s)).^2).*Se; 
Smu6_s(:,ibeta_spread)=((abs(mu6_s)).  ^2).*Se; 
Seta2_k(:,ibeta_spread)=((abs(eta2_k)).^2).*Se; 

Smu2_k(:,ibeta_spread)  =  ((abs(mu2_k)).^2).*Se; 
Seta6_k(:,ibeta_spread)=((abs(eta6_k)).^2).*Se; 

Smu6_k(:,ibeta_spread)  =  ((abs(mu6_k)).^2).*Se; 
end  %  Loop  on  beta_spreading  ends 

% 

%  Initiaiize  variabies  before  integrating 
% 

Sf_i=0; 

Sxi_s_i=0; 
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Sxi_k_i=0; 

Sxi0_s_i=0; 

Sxi0_k_i=0; 

SF2s_t_i=0; 

SF2k_t_i=0; 

Seta  2_s_i=0; 

Smu2_s_i=0; 

Seta6_s_i=0; 

Smu6_s_i=0; 

Seta2_k_i=0; 

Smu2_k_i=0; 

Seta6_k_i=0; 

Smu6_k_i=0; 

% 

%  Integrate  response  spectra  over  omega  and  theta 
% 

for  J  =2:l:beta_incr_values,  %  13  for  15  deg  increments,  37  for  5  deg  increments 
for  I  =2:l:filesize(2), 

% 

delta_omega=abs(omegae(l-l)-omegae(l)); 
deita_theta=beta_incr*  pi/180; 

% 

Sf_sum=Sf(l,J)+Sf(l-l,J)+Sf(l,J-l)+Sf(l-l,J-l); 

Sxi_s_sum=Sxi_s(l  ,J  )+Sxi_s(l-l,J  )+Sxi_s(l  ,J  -l)+Sxi_s(l-l,J  - 1); 

Sxi_k_sum=Sxi_k(l  ,J  )+Sxi_k(l-l,J  )+Sxi_k(l  ,J  -l)+Sxi_k(l-l,J-l); 
SxiO_s_sum=SxiO_s(l,J)+SxiO_s(l-l,J)+  SxiO_s(l  ,J-1)  +  SxiO_s(l-l,J-l); 
SxiO_k_sum=SxiO_k(l,J)  +  SxiO_k(l-l,J )  +  SxiO_k(l  ,J -1)  +  SxiO_k(l-l,J-l); 
SF2s_t_sum=SF2s_t(l,J)  +  SF2s_t(l- 1,J )  +  SF2s_t(l,J-l)  +  SF2s_t(l-l,J-l); 
SF2k_t_sum=SF2k_t(l,J)  +  SF2k_t(l-l,J )  +  SF2k_t(l,J-l)  +  SF2k_t(l-l,J-l); 
Seta2_s_sum=Seta2_s(l,J)  +  Seta2_s(l- 1,J )  +  Seta2_s(l,J-l)  +  Seta2_s(l-l,J-l); 
Smu2_s_sum=Smu2_s(l,J )  +  Smu2_s(l-l,J )  +  Smu2_s(l,J-l)  +  Smu2_s(l-l,J-l); 
Seta6_s_sum=Seta6_s(l,J)  +  Seta6_s(l- 1,J )  +  Seta6_s(l,J-l)  +  Seta6_s(l-l,J-l); 
Smu6_s_sum=Smu6_s(l,J )  +  Smu6_s(l-l,J )  +  Smu6_s(l,J-l)  +  Smu6_s(l-l,J-l); 
Seta2_k_sum=Seta2_k(l,J)  +  Seta2_k(l-l,J )  +  Seta2_k(l,J-l)  +  Seta2_k(l-l,J-l); 
Smu2_k_sum=Smu2_k(l  ,J )  +  Smu2_k(l-l,J )  +  Smu2_k(l,J-l)  +  Smu2_k(l-l,J-l); 
Seta6_k_sum=Seta6_k(l,J)  +  Seta6_k(l-l,J )  +  Seta6_k(l,J-l)  +  Seta6_k(l-l,J-l); 
Smu6_k_sum=Smu6_k(l  ,J )  +  Smu6_k(l-l,J )  +  Smu6_k(l,J-l)  +  Smu6_k(l-l,J-l); 
% 

Sf_i=Sf_i  +  0.25*delta_omega*delta_theta*Sf_sum; 

Sxi_s_i=Sxi_s_i  +  0.25*delta_omega*delta_theta*Sxi_s_sum; 

Sxi_k_i=Sxi_k_i  +  0.25*delta_omega*delta_theta*Sxi_k_sum; 

SxiO_s_i=SxiO_s_i  +  0.25*delta_omega*delta_theta*Sxi0_s_sum; 
SxiO_k_i=SxiO_k_i  +  0.25*delta_omega*delta_theta*Sxi0_k_sum; 
SF2s_t_i=SF2s_t_i  +  0.25*deita_omega*delta_theta*SF2s_t_sum; 
SF2k_t_i=SF2k_t_i  +  0.25*delta_omega*delta_theta*SF2k_t_sum; 
Seta2_s_i=Seta2_s_i  +  0.25*delta_omega*delta_theta*Seta2_s_sum; 
Smu2_s_i=Smu2_s_i  +  0.25*delta_omega*delta_theta*Smu2_s_sum; 
Seta6_s_i=Seta6_s_i  +  0.25*delta_omega*delta_theta*Seta6_s_sum; 
Smu6_s_i=Smu6_s_i  +  0.25*delta_omega*delta_theta*Smu6_s_sum; 
Seta2_k_i=Seta2_k_i  +  0.25*delta_omega*deita_theta*Seta2_k_sum; 

Smu2_k_l  =Smu2_k_i  +  0.25*delta_omega*delta_theta*Smu2_k_sum; 
Seta6_k_i=Seta6_k_i  +  0.25*delta_omega*delta_theta*Seta6_k_sum; 
Smu6_k_i=Smu6_k_i  +  0.25*delta_omega*delta_theta*Smu6_k_sum; 
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end 

end 

% 

%  RMS  values 
% 

RMS_f  =  sqrt(Sf_i); 

RMS_xi_s  =  sqrt(Sxi_s_i); 

RMS_xi_k  =  sqrt(Sxi_k_i); 

RMS_xiO_s  =  sqrt(SxiO_s_i); 

RMS_xiO_k  =  sqrt(SxiO_k_i); 

RMS_F2s_t  =  sqrt(SF2s_t_i); 

RMS_F2k_t  =  sqrt(SF2k_t_i); 

RMS_eta2_s=  sqrt(Seta2_s_i); 

RMS_mu2_s  =  sqrt(Smu2_s_i); 

RMS_eta6_s=  sqrt(Seta6_s_i); 

RMS_mu6_s  =  sqrt(Smu6_s_i); 

RMS_eta2_k=  sqrt(Seta2_k_i); 

RMS_mu2_k  =  sqrt(Smu2_k_i); 

RMS_eta6_k=  sqrt(Seta6_k_i); 

RMS_mu6_k  =  sqrt(Smu6_k_i); 

% 

RMS_f_vector(ibeta,iSpeed)  =  RMS_f/(rho*g*L^2); 
RMS_eta2_s_vector(ibeta,iSpeed)=  RMS_eta2_s/RMS_mu2_s; 
RMS_eta6_s_vector(ibeta,iSpeed)=  RMS_eta6_s/RMS_mu6_s; 
RMS_eta2_k_vector(ibeta,iSpeed)=  RMS_eta2_k/RMS_mu2_k; 
RMS_eta6_k_vector(ibeta,iSpeed)=  RMS_eta6_k/RMS_mu6_k; 

% 

if  RMS_eta2_s_vector(ibeta,iSpeed)  >  motion_ratio 
RMS_eta2_s_vector(ibeta,iSpeed)  =  motion_ratio; 
end 

if  RMS_eta6_s_vector(ibeta,iSpeed)  >  motion_ratio 
RMS_eta6_s_vector(ibeta,iSpeed)  =  motion_ratio; 
end 

if  RMS_eta2_k_vector(ibeta,iSpeed)  >  motion_ratio 
RMS_eta2_k_vector(ibeta,iSpeed)  =  motion_ratio; 
end 

if  RMS_eta6_k_vector(ibeta,iSpeed)  >  motion_ratio 
RMS_eta6_k_vector(ibeta,iSpeed)  =  motion_ratio; 
end 

end  %  Loop  on  beta  ends 
end  %  Loop  on  speed  ends 

% 

%  Save  results 
% 

clear  msv* 
clear  mkv* 

I_string=num2str(l); 

T_m_string=num2str(T_m); 

save_filename=strcat('h_speed_Bretch_',HS_string,'_',T_m_string,'_',l_string,'.mat'); 

save(save_filename); 
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APPENDIX  E.  MATLAB  AMPLITUDE  PLOT  ALGORITHM 


% 

%  Contour  plots  -  Bretchneider  short  crested  seas 
% 

type=input('Plot  type  (l=Speed/heading)  (2=Waveheight/heading)  =  '); 

% 

%  Speed- Heading  plots 
% 

if  type=  =  l 

HS  =input('Significant  Wave  Height  (feet)  =  '); 

T_m  =input('Modal  Period  (sec)  =  '); 

I  =input('Length  (l/L)  =  '); 

HS_string=num2str(HS); 

I_string=num2str(l); 

T_m_string=num2str(T_m); 

% 

%  Horizontal  plane 
% 

load_filename=strcat('h_speed_Bretch_',HS_string,'_',T_m_string,'_',l_string,'.mat'); 

load(load_filename); 

% 

figure(l) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,0:l:20); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold  on 

c_p=[0.000:0. 005:1.0]; 

contour(X', Y',RMS_f_  vector,  c_p),colorbar 

titleCTransverse  Connection  Force') 

% 

figure(2) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,0:l:20); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold  on 

c_p=[0.000:0.010:motion_ratio]; 

contour(X', Y',RMS_eta2_s_ vector, c_p),caxis([0  motion_ ratio] ),colorbar 
titleCLeading  Ship,  Sway') 

% 

figure(3) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,0:l:20); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold  on 

c_p=[0.000:0.010:motion_ratio]; 

contour(X',Y', RMS_eta6_s_ vector, c_p),caxis([0  motion_ ratio] ),colorbar 
titleCLeading  Ship,  Yaw') 

% 

figure(4) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,0:l:20); 

[X,Y]=pol2cart(th,r); 
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h=polar(th,r);delete(h); 
hold  on 

c_p=[0.000:0.010:motion_ratio]; 

contour(X',Y',RMS_eta2_k_vector,c_p),caxis([0  motion_ratio]),colorbar 
titleCTrailing  Ship,  Sway') 

% 

figure(5) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,0:l:20); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold  on 

c_p=[0.000:0.010:motion_ratio]; 

contour(X',Y',RMS_eta6_k_vector,c_p),caxis([0  motion_ ratio] ),colorbar 
titieCTraiiing  Ship,  Yaw') 

% 

%  Verticai  piane 
% 

ioad_fiiename=strcat('v_speed_Bretch_',HS_string,'_',T_m_string,'_',i_string,'.mat'); 

ioad(load_fiiename); 

% 

figure(6) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,0:l:20); 

[X,Y]=pol2cart(th,r); 
h=poiar(th,r);deiete(h); 
hoid  on 

c_p=[0.000:0. 005:1.0]; 

contour(X', Y',RMS_f_  vector,  c_p),coiorbar 

titieCVerticai  Connection  Force') 

% 

figure]?) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,0:l:20); 

[X,Y]=pol2cart(th,r); 
h=poiar(th,r);deiete(h); 
hold  on 

c_p=[0.000:0.010:motion_ratio]; 

contour(X',Y', RMS_eta3_s_ vector, c_p),caxis([0  motion_ ratio] ),colorbar 
titieCLeading  Ship,  Heave') 

% 

figure(8) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,0:l:20); 

[X,Y]=pol2cart(th,r); 
h=poiar(th,r);deiete(h); 
hoid  on 

c_p=[0.000:0.010:motion_ratio]; 

contour(X',Y', RMS_eta5_s_ vector, c_p),caxis([0  motion_ ratio] ),coiorbar 
titieCLeading  Ship,  Pitch') 

% 

figureO) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,0:l:20); 

[X,Y]=pol2cart(th,r); 
h=poiar(th,r);deiete(h); 
hold  on 

c_p=[0.000:0.010:motion_ratio]; 
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contour(X',Y',RMS_eta3_k_vector,c_p),caxis([0  motion_ratio]),colorbar 
titleCTrailing  Ship,  Heave') 

% 

figure(lO) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,0:l:20); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold  on 

c_p=[0.000:0.010:motion_ratio]; 

contour(X',Y',RMS_eta5_k_vector,c_p),caxis([0  motion_ ratio] ),colorbar 
titleCTrailing  Ship,  Pitch') 
end 

if  type=  =  2 

V  =input('Speed  (knots)  =  '); 

T_m  =input('Modal  Period  (sec)  =  '); 

I  =input('Length  (l/L)  =  '); 

V_string=num2str(V); 

I_string=num2str(l); 

T_m_string=num2str(T_m); 

% 

%  Horizontal  plane 
% 

load_filename=strcat('h_waveheight_Bretch_',V_string,'_',T_m_string,'_',l_string,'.mat'); 

load(load_filename); 

% 

figure(l) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,0.5:0.5:30); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold  on 

c_p=[0.000:0. 005:1.0]; 

contour(X', Y',RMS_f_  vector,  c_p),colorbar 

titleCTransverse  Connection  Force') 

% 

figure(2) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,0.5:0.5:30); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold  on 

c_p=[0.000:0.010:motion_ratio]; 

contour(X',Y', RMS_eta2_s_ vector, c_p),caxis([0  motion_ ratio] ),colorbar 
titleCLeading  Ship,  Sway') 

% 

figure(3) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,0.5:0.5:30); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold  on 

c_p=[0.000:0.010:motion_ratio]; 

contour(X',Y', RMS_eta6_s_ vector, c_p),caxis([0  motion_ ratio] ),colorbar 
titleCLeading  Ship,  Yaw') 

% 

figure(4) 
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[th,r]=meshgrid((0:beta_incr:360)*pi/180,0.5:0.5:30); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold  on 

c_p=[0.000:0.010:motion_ratio]; 

contour(X',Y',RMS_eta2_k_vector,c_p),caxis([0  motion_ratio]),colorbar 
titleCTrailing  Ship,  Sway') 

% 

figure(5) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,0.5:0.5:30); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold  on 

c_p=[0.000:0.010:motion_ratio]; 

contour(X',Y',RMS_eta6_k_vector,c_p),caxis([0  motion_ratio]),colorbar 
titleCTrailing  Ship,  Yaw') 

% 

%  Vertical  plane 
% 

load_filename=strcat('v_waveheight_Bretch_',V_string,'_',T_m_string,'_',l_string,'.mat'); 

load(load_filename); 

% 

figure(6) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,0.5:0.5:30); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold  on 

c_p=[0.000:0. 005:1.0]; 

contour(X', Y',RMS_f_  vector,  c_p),colorbar 

titleCVertical  Connection  Force') 

% 

figure]?) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,0.5:0.5:30); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold  on 

c_p=[0.000:0.010:motion_ratio]; 

contour(X',Y', RMS_eta3_s_ vector, c_p),caxis([0  motion_ ratio] ),colorbar 
titleCLeading  Ship,  Heave') 

% 

figure(8) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,0.5:0.5:30); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold  on 

c_p=[0.000:0.010:motion_ratio]; 

contour(X',Y', RMS_eta5_s_ vector, c_p),caxis([0  motion_ ratio] ),colorbar 
titleCLeading  Ship,  Pitch') 

% 

figureO) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,0.5:0.5:30); 

[X,Y]=pol2cart(th,r); 

h=polar(th,r);delete(h); 
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hold  on 

c_p=[0.000:0.010:motion_ratio]; 

contour(X',Y',RMS_eta2_k_vector,c_p),caxis([0  motion_ratio]),colorbar 
titleCTrailing  Ship,  Heave') 

% 

figure(lO) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,0.5:0.5:30); 

[X,Y]=pol2cart(th,r); 

h=polar(th,r);delete(h); 

hold  on 

c_p=[0.000:0.010:motion_ratio]; 

contour(X',Y',RMS_eta5_k_vector,c_p),caxis([0  motion_ratio]),colorbar 
titleCTrailing  Ship,  Pitch') 
end 
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